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Inclusive production of KgK$ in ep collisions was studied with the ZEUS detector. Signif- 
icant production of J PC — 2 ++ tensor mesons and of the ++ glueball candidate /o(1710) 
was seen. Masses and widths were compared with previous experiments. The HI Collabo- 
ration saw a charm pentaquark candidate in the D*p spectrum at 3.1 GeV, which was not 
confirmed by a ZEUS higher statistics search. With the full HERA statistics, HI did not 
see a signal in this region. Masses, widths and helicity parameters of excited charm and 
Q |' charm-strange mesons were measured by ZEUS. Rates of c quarks hadronising into these 

D ' mesons were determined and a search for a radially excited charm meson was performed. 
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1 Introduction 

> 

■^J" ' The HERA ep collider operated with electrons or positrons at 27.6 GeV and protons at 820 

r*"** . or 920 GeV. Each of the two general purpose experiments HI and ZEUS collected during 

1995 - 2000 ("HERA I") w 120 pb- 1 and during 2003 - 2007 ("HERA II") w 370 pb- 1 . Two 
kinematic regions have been explored: Deep inelastic scattering (DIS) with photon virtuality 



*^ | Q 2 > 1 GeV 2 , where the scattered electron is visible in the main detector and photoproduction 

(PHP) with < Q 2 3 • 10~ 4 GeV 2 , where the virtual photon radiated from the incoming 
electron is quasi-real. The sample is dominated by PHP events. 
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; 2 Glueball search in the KgKg system 

h ; 

Glueballs are predicted by QCD. The lightest glueball is expected to have J PC = ++ and a 
mass in the range 1550-1750 MeV [T] and can mix with qq scalar meson nonet 1=0 states of 
similar mass. There are four such established states: /o(980), /o(1370), /o(1500) and /o(1710), 
but only two can fit into the nonet. The /o(1710) state is considered as a possible glueball 
candidate. The KgKg system can couple to J PC = ++ and 2 ++ . Therefore, it is a good place 
to search for the lowest lying ++ glueball. 

2.1 Previous results 

The e + e~ experiments TASSO and L3 studied the exclusive reaction 77 — > KgKg. L3 [5] 
saw 3 peaks and attributed them to /2(1270)/a2(1320), f 2 {l525) and /o(1710). A maximum 
likelihood fit with 3 Breit-Wigner (BW) functions plus background yielded / 2 (1525) mass and 
width values consistent with the Particle Data Group (PDG) [T] and a 4 standard deviation 
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(s.d.) signal for /o(1710) with mass and width values above PDG. The TASSO [3] K° S K° S 
spectra had no /2(1270)/a2(1320) signal and a sizable / 2 (1525) enhancement. The result was 
interpreted by interference effects between the 3 J p = 2 + resonances /2(1270), a 2 (1320) and 
/ 2 (1525) and the spectra was fitted as a sum of 3 coherent BW functions. Based on SU(3) 
symmetry arguments [4], the sign of the a 2 (1320) term for KgKg is negative and the coefficients 
of the / 2 (1270), a 2 (1320) and / 2 (1525) BW amplitudes are +5, -3 and +2. respectively. 



2.2 This analysis 



The reaction e ± p — ► K^Kg + X was stud- 
ied [5] with the full HERA luminosity of 0.5 
fb" 1 . Both PHP and DIS events were in- 
cluded. No explicit trigger requirement was 
applied for selecting the above reaction. 

K% mesons were identified via their de- 



cay mode Kg 



A clean Kg signal 



was seen for events with > 2K S candidates. 
The number of K^Kg pairs found in the Kg 
mass range 481 < M(tt + tt~) < 515 McV is 
w 672,000. 

Figure [1] shows the K%K S mass distribu- 
tion reconstructed by combining two K s can- 
didates selected in the above mass window. 
Three peaks are seen around 1.3, 1.5 and 1.7 
GeV. No state heavier than 1.7 GeV was ob- 
served. The invariant-mass spectrum, m, was 
fitted as a sum of relativistic Breit-Wigner 
(RBW) resonances and a smoothly varying 



tj{1270)/ag(1320) 



ZEUS 0.5 tb"' 

- 4 B-W * Background 

- background 



background U (m) 



exp{—Bm), where A 



M(KX) (GeV) 

Figure 1: (a)Thc K° S K% distribution (dots). 
Solid line is the coherent fit (see text); back- 
ground function is given by the dashed line. 
(b)Background-subtracted KgK% distribution 
(dots); solid line is the fit result. 



and B are free parameters. 

Two types of fit, as performed for the re- 
action 77 -> K° S K° S by L3 g] and TASSO [3], respectively, were tried. The first fit (not 
shown) is an incoherent sum of three modified RBW resonances, R, of the form F(m) = 
-), representing the peaks / 2 (1270)/a 2 (1320), / 2 (1525) and /o(1710). Here 



M R T t 



Cr is the resonance amplitude and Mr and are the resonance mass and width, respec- 
tively. The goodness of this fit is reasonable (\ 2 /ndf = 96/95); however, the dip between the 
/ 2 (1270)/a 2 (1320) and / 2 (1525) is not well reproduced. 

Figure Q] shows a coherent fit motivated by SU(3) predictions |3]. Each resonance amplitude, 



R, is described by the RBW form [5] BW{R) 



M n 



The decays of the tensor 



M^-m 2 -iM R r R ■ 

2 + ) mesons / 2 (1270), a 2 (1320) and / 2 (1525) into the two pseudoscalar ( J p = + ) mesons 
K°K° are related by SU(3) symmetry with a specific interference pattern. The intensity is the 
modulus-squared of the sum of these 3 amplitudes plus the incoherent addition of /o(1710) and 
a non-resonant background. 

Assuming SU(3) symmetry and a direct coupling of the 2 + states to the exchanged photon, 
the fitted function to the m(K^K^) spectra is given by F(m) = o[5 • BW(f 2 {1270)) - 3 ■ 
£1^(02(1320)) + 2 • BVF(/ 2 (1525))] 2 + 6[ J BH / (/ (1710))] 2 + c ■ U{m), where a,b,c as well as the 
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resonance masses and widths were free parameters in the fit. The background-subtracted mass 
spectrum is shown in FigQJb). The fit quality is good (x 2 /ndf = 86/97). The peak around 1.3 
GeV is suppressed due to the destructive interference between /2(1270) and a 2 (1320) and the 
dip between / 2 (1270)/a 2 (1320) and / 2 (1525) is well reproduced. The number of fitted / (1710) 
events is 4058 ± 820, which has « 5 s.d. significance. Its mass is consistent with a J ° = ++ 
glucball candidate, but it cannot be a pure glueball if it is the same state as in 77 — - > K^Kg. 



Fit 


No interference 


Interference 


PDG 2007 Values 


X z /ndf 


96/95 


86/97 


in McV 


Mass 


Width 


Mass 


Width 


Mass 


Width 


/ 2 (1270) 


1304 ± 6 


61 ± 11 


1268 ± 10 


176 ±17 


1275.4 ±1.1 


185.2+^ 


a 2 '(1320) 


1257 ±9 


114 ±14 


1318.3 ±0.6 


107 ±5 


#(1525) 


1523 ±3i^ 


71 ±5±£' 


1512 ± 3^ 6 


83 ±9tl 


1525 ±5 


73t* 


/o(1710) 


1692 ± 6±% 


125 ± 12+^ 


1701 ±5+3 


100 ± 24^ 9 


1724 ±7 


137 ±8 



Table 1: Fitted masses and widths for / 2 (1270), a 2 (1320), / 2 (1525) and /o(1710) from the inco- 
herent and coherent fits compared to PDG. The first error is statistical. For / 2 (1525), /o(1710) 
the second errors are systematic uncertainties. 
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The masses and widths obtained from 
both fits are shown in Table Q] and com- 
pared to PDG pp. The no-interference 
fit yields a narrow width for the combined 
/ 2 (1270)/a 2 (1320) peak, as was also seen by 
L3 [5] . The fit with interference yields widths 
close to the PDG values for all observed reso- 
nances. The a 2 (1320) mass is below the PDG 
value. The / 2 (1525) and /o(1710) masses 
are somewhat below PDG with uncertain- 
ties comparable with the PDG ones. A fit 
without /o(1710) is strongly disfavoured with 
X 2 /ndf = 162/97. 

3 Charm pentaquark search 
in the D*p system 

A narrow exotic baryon with strangeness +1 
around 1530 MeV decaying into K + n or 
K p was seen by various experiments and 
attributed to the G + = uudds pentaquark 
state predicted by Diakonov et al.[B]. If 
a strange pentaquark exists, charmed pen- 

taquarks, 9° = uuddc, could also exist. If M(6°) > M(D*) + M{p) = 2948 MeV, it can 
decay to D* ± p T . 

The HI Collaboration saw [7] in a DIS HERA I sample of f 
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Figure 2: M(D p^) from HI DIS HERA I, 
compared with fit results where both signal and 
background components are included and where 
only background is included. 
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a narrow signal of 50.6 ± 11.2 events in the D* ± p T invariant mass at 3.1 GeV (FigJ^J) with a 
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width consistent with the mass resolution and a rate of w 1% of the visible D* production. 

ZEUS searched for a 0° signal in the D^p* mode with the full HERA I PHP + DIS 
data sample (SJ. Clean D*^ signals were seen in the AM = M(L>* ± ) - M(D°) plots. Two 



D 



*± 



ZEUS 



D°TTg decay channels were used with D° ->■ i^ =F 7r ± and D° -> A" T 7r ± 7r + 7r". The 0° 
search was performed in the kinematic range |7y(D*)| < 1.6 and pr(D*) > 1.35(2.8) GeV and 
with AM values between 0.144 - 0.147 (0.1445 - 0.1465) GeV for the Kirir (Kiriririr) channel. 
In these bands sa 62000 D*'s were obtained after subtracting wrong-charge combinations with 
charge ±2 for the D° candidate and ±1 for the D* candidate. Selecting DIS events with 
Q 2 > 1 GeV 2 yielded smaller, but cleaner D* signals with w 13500 £>*'s. 

Protons were selected with momentum 
P(p) > 0.15 GeV. To reduce the pion and 
kaon background, a parameterisation of the 
expected dE/dx as a function of P/m was 
obtained using tagged protons from A decays 
and tagged pions from K s decays. The x 2 
probability of the proton hypothesis was re- 
quired to be above 0.15. 

Figure [3] shows the M(D* ± p T ) distribu- 
tions for the D° -> Kir (left) and D° -> 
Kiririr (right) channels for the full (up) and 
the DIS (down) samples. No narrow signal is 
seen in any of the distributions. 95% C.L. up- 
per limits on the fraction of D* mesons origi- 
nating from 0° decays, R(Q° -► D*p/D*), 
were calculated in a signal window 3.07 < 
M(D*p) < 3.13 GeV for the Kmr and 
Ktttttttt channels. The M(D*p) distributions 
were fitted to the form x a e~ bx+cx , where x = 
M{D*p) - M (D*) - m p (PDG). The number 
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Figure 3: 



MiD^pT) 



from ZEUS HERA I. 
Solid curves are fits to a background function. 
Shaded historgams are MC 0° signals, nor- 
malised to 0°/Z?* = 1%, on top of the back- 
ground fit. 



of reconstructed 6° baryons was estimated 
by subtracting the background function from 
the observed number of events in the signal 
window, yielding i?(9° -► D*p/D*) < 0.23% 
and < 0.35% for the full and DIS combined 
two channels. A visible rate of 1% for this 

fraction is excluded by 9 s.d. (5 s.d.) for the full (DIS) combined sample. The acceptance- 
corrected rates are, respectively, 0.37% and 0.51%. The 95% C.L. upper limit on the fraction 
of charm quarks fragmenting to 0° times the branching ratio 0" — > D*p for the combined two 
channels is f(c -► 0°) • B e o_ D . p < 0.16% (< 0.19%) for the full (DIS) sample. 

In a HERA II DIS data sample that is f=a 4 times larger than the HERA I sample, HI does 
not see any significant peak at 3.1 GeV (Figj4j) . A preliminary 95% C.L. for the ratio of D*p 
to D* is 0.1%. 



4 Excited charm and charm- strange mesons 

The large charm production at HERA allows to search for excited charm states. ZEUS studied 
the orbitally excited states L>i(2420)° -> D* ± 7r =F (J p = 1+), L>^(2460)° -► D* ± tt^,D ± tt^ 
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D^K^D^K* {J 1 



1 + ) and searched for the radially 



(J p = 2+) and D s i(2536) ± 

excited state L>*'(2640) ± -> D^n+n' (J p = 1" ?) with a HERA I PHP + DIS sample[9]. 

A large sample of events has been col- 
lected with the ground state charm mesons 
D* ± ,D a ,D ± . The number of D* ± mesons 
was obtained by subtracting the wrong 
charge background. The number of — > 
A'T 7r ± 7 r ± and D°{D°) K^ir* was ex- 
tracted from fits to a modified Gauss func- 
tion, Gauss mod ~ exp(~0.5x 1+ ( 1 + a - 5 ^ ), 
where ir = (M — Mu) /<r, plus a back- 
ground function. For the D* , both D° decay 
modes to AT7T and Kmnr were used. 

4.1 Excited charm mesons 

To reconstruct the excited charm mesons, a 
D*^ or Z) ± candidate was combined with a 




M(D*p) [ GeV ] 



Figure 4: M{D* ± p T ) from HI DIS HERA II. 
The solid line is a background parametrisation. 
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pion of opposite charge, n a . Figure [5] shows 
the "extended" mass difference distributions 
M{D* ± -K a )- M(D* ± )+M{D*) PDG (upper plot) and M (D ± ir a ) 
plot). A clear excess is seen in M(D* ± ttJ) around the D^/Df 
near the D^ mass is seen in M{D 7rJ). No excess is seen for wrong charge combinations, 
where D*(D) and 7r a have the same charge. 

To distinguish between the Dj and 
-D2 , the helicity angular distribution, 
parametrised as dN/d cos a ~ 1 + h cos 2 a, 
was used. Here a is the angle between the ir a 
and 7Ts momenta in the D* rest frame. The 
helicity parameter h is predicted [TU] to be 
3(-l) for pure D-wave D? (£>|°). 

Figure [5] shows the D^ita "extended" 
mass difference in 4 helicity | cosa| intervals. 
The D® contribution increases with |cosa| 
and dominates for |cosa| > 0.75. A si- 
multaneous fit was performed to the 4 helic- 
ity regions of Fig[6] and to the M(Dtt) dis- 
tribution of FigEl The data is described 
well with 15 free parameters (signal yields, 
masses, D® width and helicity). The fitted 
masses agree with PDG. The fitted D® width 
is 53.2±7.2(stat.)±lll(syst.) McV compared 
to 20.4 ± 1.7 MeV of PDG. The fitted D? he- 
licity (5.9 + ?'7(stat.) +2 'n(swst.)) is consistent , > , , > 

v '~ 1 - oy y " Figure 5: M B irj and M(D ± n a ) distribu- 

with a pure D-wave. . v . . ' 

tions. Solid curves arc simultaneous fit; dashed 

curves are background; histograms are wrong- 
charge combinations. 
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4.2 Excited charm strange mesons 



To reconstruct the — ► D* ± K 1 ^ decays, 
a candidate was formed by combining 
a D* candidate with a reconstructed Kg of 

the same event. Figure [7] (upper plot) shows the "extended" mass difference distribution 
M(D* ± A'g) - M(D*±) + M(D*) PDG + M(K°) PDG . A clear D s i(2536) ± signal is seen. 
The decay mode Df± — > D*°K^ is reconstructed from the "extended" mass difference 
M(D°K a ) - M(D°) + AI(D") PDG . A nice Df x signal is seen (Figure lower plot) at a mass 
shifted down by w 142 MeV from the mass. The signal is a feed-down from 
with Z?*° — > D°ir°, D 0r y. An unbinned likelihood fit was performed using simultaneously values 
of M(D°K a ), M{D* ± K%) and cos a for the D r± K^ combinations. Yields and widths of both 
signals and the mass and helicity parameter were free parameters of the fit. The fitted 
D s i helicity parameter is fo(-D^) = —Q.7A'^^{stat.)^ t ^{syst.). It is inconsistent with a pure 
J p — 1 + D-wave and is barely consistent with a pure J p — 1 + S-wave, indicating a significant 
S — D mixing. 

The helicity angular distribution form of a 1 + state for any D- and S-wave mixing is: 
dN/dcosata r+(l-r)(l+3 cos 2 a)/2+^2r(l - r) cos 0(1—3 cos 2 a), wherer = rs/(r s + T D ), 
Tg/u is the S/D wave partial width and </> is relative phase between the 2 amplitudes, cos0 = 
^ 3 ?^/2 ^^S ure El shows a range, restricted by the measured ft(D^) value and its un- 

certainties, in a plot of cos0 versus r. The measurement suggests a significant contribution of 
both D- and S-wave amplitudes to the D sl (2536) ± -> D*±Kl decay. The ZEUS range agrees 
nicely with the BELLE result and roughly with the CLEO measurement. 
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Figure 7: M( J D* ± X°) and M(D°K ± ) distribu- 
tions. Solid curves are simultaneous fit; dashed 
curves are background. 
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Figure 6: M(D 7r a ) distributions in 4 helicity 
intervals. 
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5 Branching ratios and fragmentation fractions 



Using the ZEUS measured fractions /(c — ► D* + ) 
and f(c — > I?" 1 ")!!!], the following decay 

rate ratios were derived: „ — - = 

2.8± 0.8±°^ (PDG: 2.3±0.6); D t^ D '° K+ = 

D+-»D* + K° 

2.3 ± 0.6 ± 0.3 (PDG: 1.27 ±0.21). 

Assuming isospin conservation for D® and 

D 2° and B D+_->D*+K° + B D+ 1 ^D* a K+ = 1 

yields a strangeness suppression of excited D 
mesons /(c -► D+)/f{c -► D?) = 0.31 ± 
0.06(stat.)t°;^( S2/s i.). 

In Table [5] the ZEUS fragmentation frac- 
tions of the excited charm mesons are com- 
pared with e+e" values. The results are con- Figure g . CQS(j) yg _ r s /(r s + L D ) for 
sistent within errors. D + ^ D *+ R o decay _ 

DELPHI saw a narrow peak in D* ± tt + tt 
at 2637 MeV [12] and attributed it to a ra- 
dially excited D* ± . No signal was seen in 
ZEUS and a 95% C.L. upper limit of /(c -> 

D*^) • B D «< ± ^ D , +7r+7r - < 0.4% was set, compared to the weaker limit of OPAL (0.9%) [13]. 





f(.c->Dt)[%] 


f(c^D?)[%] 


/(c 


ZEUS 


3.5±0.4i^ 


3.8±0.7i^ 


1.11 ± 0.16t£¥o 


OPAL 


2.1 ±0.8 


5.2 ±2.6 


1.6 ±0.4 ±0.3 


ALEPH 






0.94 ±0.22 ±0.07 



Table 2: The fractions of c quarks hadronising into D®, D^ and Dji mesons. 
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